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ABSTRACT 

Corrections to the magnitudes of high redshift objects due to intergalactic attenuation 
are computed using current estimates of the properties of the intergalactic medium. 
The results of numerical simulations are used to estimate the contributions to resonant 
scattering from the higher order Lyman transitions. Differences of 0.5 — 1 magnitude 
from the previous estimate of Madau (1995) are found. Intergalactic fcicM-corrections 
and colours are provided for high redshift starburst galaxies and Type I and Type II 
QSOs for several filter systems used in current and planned deep optical and infra-red 
surveys. 
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1 INTRODUCTION 

Over the past decade, deep optical and infra-red surveys 
have enabled giant strides to be taken in elucidating the na¬ 
ture and properties of objects that populate the high redshift 
universe. The band-dropout method has unveiled a popu¬ 
lation of Lyman break galaxies at « « 3 (Guhathakurta 
et al. 1990; Bithell 1991; Steidel & Hamilton 1992, 1993; 
Steidel, Pettini & Hamilton 1995; Steidel et al. 2003) and 
higher (Sawicki & Thompson 2005). The selection method 
was successfully applied to the Hubble Deep Field (HDF) 
(Giavalisco, Steidel & Macchetto 1996; Steidel et al. 1996), 
broadening the redshift range and volume coverage over pre¬ 
vious surveys. Most recently the Ultra Deep Field (UDF) 
was exploited to discover objects as distant as 2 6 (Stan¬ 

way, Bunker & McMahon 2003; Bouwens et al. 2004; Yan 
& Windhorst 2004; Giavalisco et al. 2004). 

An alternative selection method for identifying high 
redshift objects relies on combinations of broadband colours 
to estimate photometric redshifts, for which the most likely 
redshift is assigned based on predicted spectral energy dis¬ 
tributions (Sawicki, Lin & Yee 1997; Gsabai et al. 2000). 

The modelling of the high redshift objects through 
population synthesis, applied to a combination of spectro¬ 
scopic data and broad-band colours, suggests that most of 
the high redshift objects are star-forming galaxies (Madau 
et al. 1996; Metcalfe et al. 2001; Papovich, Dickinson & Fer¬ 
guson 2001; Pettini et al. 2001; Shapley et al. 2003). A few 
of the Lyman break galaxies contain Active Galactic Nuclei 
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(Steidel et al. 2002), used to determine the faint end of the 
QSO luminosity function at 2 « 3 (Hunt et al. 2004). 

Parallel to these surveys have been several searches for 
high redshift Quasi-Stellar Objects (QSOs). The discovery 
of a few dozen 2 > 3.6 QSOs by the Sloan Digital Sky Survey 
(SDSS) has made possible a new evaluation of the bright end 
of the QSO luminosity function and its evolution at these 
high redshifts (Fan et al. 2001). The results are currently 
being revised (Richards et al. 2005) based on the much 
larger numbers now detected, including over 500 at 2 > 4 
(Schneider et al. 2005). 

Similar surveys are expected to continue well into the 
future, including some now in progress, such as the Ganada- 
France-Hawaii Telescope Legacy Survey (GFHTLS^) and 
the UKIRT Infrared Deep Sky Survey (UKIDSsk Hewett 
et al. , in preparation), or planned for telescopes currently 
under development, such as the Visible & Infrared Survey 
Telescope for Astronomy (VISTA§) and the James Webb 
Space Telescope (JWST^). 

Crucial to all these analyses is an accurate estimate of 
the amount of intergalactic attenuation due to intervening 
absorption systems. Most have relied on the standard work 
of Madau (1995), whose assessment was based on the then 
most current understanding of the properties and distribu¬ 
tions of intervening systems. Madau (1995) estimated the 
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blanketing due to the resonant (Lyman series) scattering of 
photons assuming idealised forms for the H i column density 
distribution of the absorbers. The contributions of the differ¬ 
ent orders in the Lyman series depend on the full line-shape 
of the absorber, precluding a direct scaling of line-centre 
optical depths based on pure atomic physics considerations. 
Instead the broadening of the absorbers must be included, 
the distribution of which has since been shown to be sensi¬ 
tive to column density (Kirkman & Tytler 1997; Kim, Cris- 
tiani & D’Odorico 2002b). Madau (1995) adopted a constant 
Doppler parameter for all absorption systems, varying the 
value to probe the sensitivity of the total amount of attenua¬ 
tion to this variable. A blanketing formalism based on Pois¬ 
son placement of the absorbers was used to predict the effec¬ 
tive optical depths, although in principle small scale cluster¬ 
ing of the absorption systems will affect the total amount of 
blanketing, and such correlations have been detected (Kirk¬ 
man & Tytler 1997; Kim et al. 2002a). 

Since Madau’s seminal work, numerical simulations of 
the intergalactic medium have yielded results matching 
the measured distributions of the Lya flux distributions 
to an accuracy of a few percent, as well as the line pa¬ 
rameters (allowing for extra heating) (Meiksin, Bryan & 
Machacek 2001). Numerical simulations have also repro¬ 
duced the H I column density dependence of the Doppler 
parameter envelope (Misawa et al. 2004), as well as cor¬ 
relations in the H i flux distribution (Croft et al. 2002; 
Meiksin & White 2004). The understanding of the mean 
Lya intergalactic optical depths has improved substantially 
over the past decade (see Meiksin & White 2004 for a sum¬ 
mary and Kirkman et al. 2005 for subsequent results). The 
simulations contain the information necessary to extract the 
contributions from all higher order transitions to the blan¬ 
keting. Although simulations do not recover the full num¬ 
bers of Lyman Limit Systems observed (Gardner et al. 1997; 
Meiksin & White 2004), an assessment of their numbers over 
the redshift range 0.3 ^ z 4 has been made by Stengler- 
Larrea et al. (1995). These improvements now permit a 
much more secure determination of the amount of inter¬ 
galactic attenuation to be made. 


2 NUMERICAL SIMULATION PREDICTIONS 

FOR INTERGALACTIC ATTENUATION 

Attenuation due to intervening intergalactic hydrogen arises 
through two principle mechanisms, resonant scattering by 
Lyman transitions and photoelectric absorption. Additional 
contributions are made by intervening metal systems and in¬ 
tervening helium. The contributions of metals and He i are 
small (Madau 1995), while He ii will contribute only at 
wavelengths A < 228(1 -I- z)k. for a source at redshift z. 

The hydrogen attenuation values are based on the sim¬ 
ulation of a ACDM model by Meiksin & White (2004) with 
Dm = 0.3, Q.V = 0.7, h = 0.7, n = 0.95, and normalised to 
(T8 = 0.92, consistent with WMAP constraints. The simu¬ 
lation used a pure particle mesh (PM) scheme (Meiksin & 
White 2001), mimicking the temperature of the gas using a 
polytropic equation of state and assuming the gas and dark 
matter have the same spatial distribution. The mean trans¬ 
mitted Lya fluxes are normalised using the values listed in 
Table 3 of Meiksin & White (2004). The attenuation values 


Table 1. Ratio of Tn/ro for Lyman transition n —> 1 for n = 3 
to 9. Each transition contributes at wavelength A for an object at 
redshift z only under the condition Zn < 2 , where Zn = (A/A„) —1 
and Xn is the (rest) wavelength of the transition. 


n Tn/rc 


3 

0.348(0.25(1 -1- 2 , 


< 3); 0.348(0.25(1 -1- 

{Zn 

>3) 

4 

0.179(0.25(1 -1- 2 , 

(2„ 

< 3); 0.179(0.25(1 -1- 

Zn)y-^^, (Zn 

>3) 

5 

0.109(0.25(1 -1- 2 , 

n)]D3, (2„ 

< 3); 0.109(0.25(1 -1- 

2„)(3/3, (2„ 

>3) 

6 

0.0722(0.25(1 + , 

2„))l/3 




7 

0.0508(0.25(1 + , 

2„))l/3 




8 

0.0373(0.25(1 + , 





9 

0.0283(0.25(1 + , 

2„))l/3 





resulting from higher order Lyman transitions are computed 
directly from the simulation results. 

The mean optical depth value for the Lyman transition 
n —> 1 is dehned by 

r„ = - ln(exp(-r„)), (1) 

where (exp(—r„)) is the corresponding mean transmitted 
flux. (Here, the convention Ta = T 2 , t/b = ts, etc. will be 
adopted.) For 2 < 2 < 4, the measured values of the mean 
Lya transmitted flux are fit to within 2% using 

Tc = 0.00211(1 -h 2 )® '^ iz<4). (2) 

The sharp reduction in flux at higher redshifts is adequately 
fit by 

Tc = 0.00058(1 -h zf-^ (2 > 4). (3) 


For 2 > 6, the measured values become quite uncertain, but 
the objects become dimmed by 4 magnitudes or more so 
that accurate values will normally not be required. 

The higher order terms in the Lyman series are char¬ 
acterised by the ratio Tn/ra to factor out most of the red¬ 
shift dependence. A residual dependence, however, remains. 
These are accurately fit (to within a few percent) by weak 
powers of (1 -f 2 ). The ratios for the first 7 orders (Ly/3 to 
Ly^ ) after Lya are provided in Table 1. Higher orders 
(n > 9) are found to scale very nearly according to the 
atomic physics prediction for line-centre optical depth val¬ 
ues: 


^ _ 720 

fe n{rp — 1) 


( 4 ) 


Terms up to n = 31 are included: higher order Lyman series 
transitions contribute negligibly. 

The contribution due to photoelectric absorption is split 
into two parts, the contribution from systems optically thin 
at the Lyman edge and the contribution from Lyman Limit 
Systems. The contribution from optically thin systems is 
given in the linear approximation tl = r/ro (Zuo 1992), 
where r is the proper distance between the emitting object 
at redshift 2 and the redshift zl = A/An — 1, where \l 
is the wavelength at the Lyman edge (An = 912A), A is 
the observed wavelength, and ro is the attenuation length 
of photons at the photoelectric edge. Using the result for 
the attenuation length from Meiksin & White (2004), this is 
well-approximated by 
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Figure 1. Intergalactic transmission as a function of observed 
wavelength, (a) Mean intergalactic transmission for a source at 
redshift z = 3.5 (solid lines). The mean intergalactic transmission 
assuming no Lyman Limit Systems lie along the line-of-sight is 
shown by the dashed lines, (b) Mean intergalactic transmission 
for, from left to right, sources at redshifts z = 2.5, 3.5 and 4.5. 
The heavy lines are the estimates from this paper. In both panels, 
the light lines are the estimates from Madau (1995). 


contribute only a small amount while very optically thick 
systems are too few in number to contribute much on av¬ 
erage. While the photoelectric contribution due to the opti¬ 
cally thin IGM is well-determined due to the large num¬ 
ber of optically thin systems, the contribution from Ly¬ 
man Limit Systems is highly variable, depending on the 
chance that a Lyman Limit System lies along the relevant 
path length. For instance, for a source a,t z = 3.5, the 
Lyman limit lies at the observed wavelength A = 4104A. 
The number of Lyman Limit Systems that will fall at a 
redshift corresponding to the [/„-band wavelength range of 
approximately SOOOA—4000A is 2.1. The Poisson probabil¬ 
ity that no Lyman Limit System lies in this range is then 
= 0.12. This will introduce large fluctuations in the 
amount of attenuation at these short wavelengths (Zuo & 
Phinney 1993; Madau 1995), a fact which should be borne 
in mind when interpreting the number of source detections 
in Lyman dropout surveys. 

An indication of the level of contribution of Lyman 
Limit Sytems to the total optical depth is provided in Fig¬ 
ure la, which shows the level of intergalactic transmission as 
a function of observed wavelength and the transmission with 
the contribution of Lyman Limit Systems removed. The to¬ 
tal intergalactic transmissions allowing for H i attenuation 
are shown as a function of the observed wavelength for 
sources at z = 2.5, 3.5 and 4.5 in Figure lb, and in Table 2. 
The figures include comparisons with the corresponding es¬ 
timates of Madau (1995), which tend to lower transmission 
levels, primarily as a result of differences in the estimates of 
the contributions of resonant absorption. 


+«)’(T 4 ^-rb)- ® 

The contribution due to Poisson-distributed Lyman 
Limit Systems is given by (Zuo 1992) 


LLS 


tl 



d'^N 

OtlOz' 



where d'^N/O tlOz = + zC is assumed, correspond¬ 

ing to a number density dN/dz = No{l -I- zy for systems 
with tl > 1, where Nq = A/{fi — 1). The integral has the 
power series solution 


LLS 

Tl 


No 


4+7-3/3 
X 


r(2-/l,i)-e-i-E-o4^^^ 

-3(/3-1)+7-)-1 t A \3(/3-l) / A 


,1+ 

_ ^ y'oo _ f3 — l _ ( —1)^ 

l\0Zjn=i n! 


( 7 ) 


The first 10 terms of each series provide a high level of con¬ 
vergence. The results are normalised to Nq = 0.25, /3 = 1.5 
and 7 = 1.5 (Stengler-Larrea et al. 1995). Although the 
result of Stengler-Larrea et al. applies only for 2^/4, the 
redshift dependence is extrapolated to higher redshifts here, 
noting that the contribution of these systems is uncertain at 
these redshifts. 

The total optical depth is the sum of the resonant and 
photoelectric contributions. The photoelectric contribution 
is dominated by sytems with « 1; optically thin sytems 


3 EFFECT OF INTERGALACTIC 

ATTENUATION ON BROADBAND 

MAGNITUDES 

In the AB-magnitude system, the apparent magnitude of a 
source of intrinsic flux fi, (in cgs units) measured through a 
Tilter with (normalised) transmissivity T{u) is given by 

TTiAs =-2.51ogiQ J exp(-reff)r(iA) - 48.59. (8) 

The difference between the magnitudes with and without in¬ 
tergalactic attenuation will be designated as the intergalac¬ 
tic k-correction fcicM = mAs(Teff) — mAB(reft = 0). Because 
source spectra are typically slowly varying over a bandwidth, 
the fciGM-correction is fairly independent of the nature of the 
source. The exception is when either the source has narrow 
features that dominate the light, such as strong emission or 
absorption lines, or when the amount of intergalactic atten¬ 
uation varies rapidly within a band, reducing the wavelength 
range of the source spectrum that contributes to the total 
magnitude, and so enhancing the effect of any differences 
between source spectra. 

To indicate the typical role intergalactic attenuation 
plays on the colours of high redshift objects, the effect of in¬ 
tergalactic attenuation is evaluated for four model sources: 
starbursts of ages 3 Myr and 600 Myr, roughly brack¬ 
eting the range inferred for Lyman Break Galaxies (Pa- 
povich, Dickinson & Ferguson 2001), and Type I and Type 
II QSOs. The starburst spectra assume continuous star- 
formation with solar metallicity and a Salpeter Initial Mass 
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Table 2. Intergalactic H I transmission as a function of observed wavelength. 


A (A) 

2 = 1.5 


2 = 2.5 


2 = 3.5 


2 = 4.5 


2 = 5.5 

q 

II 

q 

II 

U! 

0 

II 

U! 

1730.0 

0.323241 

0.150055 

0.074558 

0.038532 

0.020407 

0.010983 

0.005976 

0.003277 

0.001807 

0.001001 

0.000556 

0.000310 

1731.0 

0.323538 

0.150066 

0.074510 

0.038482 

0.020369 

0.010956 

0.005958 

0.003266 

0.001800 

0.000996 

0.000553 

0.000308 

1732.0 

0.323836 

0.150077 

0.074462 

0.038433 

0.020330 

0.010929 

0.005940 

0.003254 

0.001793 

0.000992 

0.000550 

0.000306 

1733.0 

0.324135 

0.150089 

0.074414 

0.038383 

0.020292 

0.010902 

0.005922 

0.003243 

0.001785 

0.000987 

0.000548 

0.000304 

1734.0 

0.324436 

0.150101 

0.074367 

0.038334 

0.020253 

0.010875 

0.005905 

0.003231 

0.001778 

0.000983 

0.000545 

0.000303 


Note: The full table is published in the electronic version of the paper. A portion is shown here only for guidance regarding its form and content. 



redshift 


Figure 2. Intergalactic attenuation fcioM-correction and colours 
for a 600 Myr starburst with continuous star formation, (a) From 
left to right, feiGM for Un, G, TZ, I and z' bands (solid lines). For 
the Un and G bands, fciGM assuming the absence of Lyman Limit 
Systems is shown by the dashed lines. The absence of Lyman 
Limit Systems leaves the results for the remaining bands unal¬ 
tered. The light solid lines correspond to the predictions for fciGM 
using the attenuation model of Madau (1995). (b) The Un — G 
(solid) and G — TZ (long dashed) colours of the starburst. The dot¬ 
ted short-dashed lines and dotted long-dashed line show Un — G 
and G — TZ, respectively, assuming the absence of Lyman Limit 
Systems. The light lines correspond to the predictions based on 
the attenuation model of Madau (1995). 


Function and were generated by the STARBURST99 model 
of Leitherer et al. (1999). The Type I QSO spectrum is the 
composite spectrum constructed from over 2200 spectra ho¬ 
mogeneously selected from the SDSS QSO survey covering 
the redshift interval 0.044 < z < 4.789 and restframe wave¬ 
length range 800A—8555A (Vanden Berk et al. 2001). The 
Type II QSO spectrum is based on observations of CXO 52 
(Stern et al. 2002), as described in Meiksin (2005). This 
Type II QSO was chosen because it has high equivalent 
width emission lines that will dominate the magnitude of 
the band in which they lie, and so represents a class of ob¬ 
jects that will have unusual broadband colours that vary 
substantially with redshift. 


The effect of intergalactic attenuation on broadband 
magnitudes is computed for several filter systems relevant 
to recent and planned major surveys: the Sloan u'g'r'i'z' 
system (assuming airmass 1.3 response for a point source); 
the UnGTZ system of Steidel’s group, supplemented by the 
7-band adopting the Harris I filter; the Hubble Space Tele¬ 
scope (HST) ACS/NICMOS B 435 (T435IT), Vme {F606W), 
1775 iF775W), 2850 {F850LP), Jno (F’llOIT) and Hieo 
{Fl&OW) bands (including instrumental responses); the 
UKIDSS ZYJHK filters; and the JWST F070W (desig¬ 
nated mo 7 o and approximated as a square transmission re¬ 
sponse) and T150IT (designated miso) bands. All magni¬ 
tudes are computed on the AB system. 

The fciGM-corrections and colours for the Un, G and TZ 
filters of Steidel & Hamilton (1992) as well as the I and 
z' bands, are shown in Figure 2 for the 600 Myr old star- 
burst. The effect of internal redenning is not included. This 
will negligibly affect fciGM, so the magnitude corrections to 
the colours due to redenning will simply be additive. Also 
shown are the results assuming no Lyman Limit Systems 
intercept the line-of-sight. The mean contribution of Lyman 
Limit Systems has a substantial effect on t/„, a small effect 
on G, and a negligible effect on the longer wavelength bands 
for fciGM < 3. A comparison with the attenuation model of 
Madau (1995) shows a reduction in fciGM of about 0.5 magni¬ 
tude, although the effect on the colours approach a shift of a 
full magnitude. In terms of photometric redshifts, the differ¬ 
ence in attenuation models produces a shift of A 2 = 0.1 —0.2 
for Un — G and G — TZ. The fciGM-corrections and various 
colour combinations for all the sources considered are shown 
in Figures 3 and 4. The fciGM-corrections are nearly indepen¬ 
dent of the source spectrum, except for the Type H QSO for 
which the magnitudes may be dominated by emission lines 
(Meiksin 2005). The colours also are similar except for the 
Type H QSO, although certain colours, such as UKIDSS 
J — K and H — K, are particularly effective at separating 
the other sources. 


4 SUMMARY 

The amount of the attenuation of light from high redshift ob¬ 
jects on passing through the IGM has been computed based 
on recent measurements of the mean transmitted Lya flux 
through the IGM and recent assessments of the numbers of 
Lyman Limit Systems. The properties of the IGM required 
to compute the attenuation due to resonant Lyman photon 
scattering were based on numerical simulations matching the 
measured properties of the IGM. Differences from the pre¬ 
dictions of the model of Madau (1995) are found for Isigm 
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Figure 3. Intergalactic attenuation feiGM-coirections as a function of source redshift for a 600 Myr starburst (solid 
lines), a 3 Myr starburst (long-dashed lines), a Type I QSO (dot-dashed lines) and a Type II QSO with high 
equivalent width emission lines (dotted lines). There is little variation between the sources except for the Type II 
QSO because of the presence of strong emission lines. The rows correspond to the following filter systems. First 
row: Sloan; second row: Steidel UnGTZ filters and /-band; third row: //5T bands i? 435 , Vgoei *775 JWST mo 7 o; 
fourth row: Sloan 2 ', HST 2350 and Jno and UKIDSS Y. The fciGM-corrections are negligible for the HST 
JWST mi 5 o and UKIDSS J, H and K bands for 2 < 8.5. 


and colours of 0.5 — 1 magnitude for bands containing rest- 
frame Lya and shorter wavelengths. 

Intergalactic fciGM-corrections are provided for filter 
systems used in current or planned deep optical and infra¬ 
red surveys, viz. the Sloan u'g'r'i'z' system, Steidel UnGTZ 
and the /-band, UKIDSS ZYJHK, the HST B 435 , Hoe, 
* 775 , 2850 , ,/iio and i/ieo bands, and the JWST FQ7QW and 
F150W bands. 

Colours based on the above bands are provided for star- 
bursts of ages 3 Myr and 600 Myr, typical of 2 « 3 Lyman 
Break Galaxies, and Type I and Type II QSOs, over the red¬ 
shift range 2.5 < 2 < 8.5. The results show which colours 
are most effective for distinguishing between different ob¬ 


jects and the typical values to expect, which may be used 
for the planning and analysis of current and upcoming deep 
surveys. 
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Figure 4. Colour predictions including the effects of intergalactic attenuation for a 600 Myr starburst (solid lines), a 3 Myr starburst 
(long-dashed lines), a Type I QSO (dot-dashed lines) and a Type II QSO with high equivalent width emission lines (dotted lines). In the 
third and fifth rows, the B, V, i, z, J and H magnitudes refer to the HST S 435 , Vgoe? ^ 775 ? ^850? '^110 -^^160 bands, respectively. The 

fourth row magnitudes refer to the UKIDSS ZYJHK bands. 
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